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Abstract 



CD ■ 

^ ' Using Tyablikov's decoupling approximation, we calculate the initial suppression rate of the Neel temperature, _Ris = — lim^^o T N 1 

• , dTyt/dx, in a quasi two-dimensional diluted Heisenberg antiferromagnet with nonmagnetic impurities of concentration x. In order 

• ' to explain an experimental fact that -Rig"' = 3.4 of the Zn-substitution is different from -Rig 18 ' = 3.0 of the Mg-substitution, 
, we propose a model in which impurity substitution reduces the intra-plane exchange couplings surrounding impurities, as well 

(~j as dilution of spin systems. The decrease of 12% in exchange coupling constants by Zn substitution and decrease of 6% by Mg 
I substitution explain those two experimental results, when an appropriate value of the interplane coupling is used. 
*T3 '. © 2008 Elsevier B.V. All rights reserved. 
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The discovery of high-T c compounds renewed interests 
in quasi two-dimensional quantum antiferromagnets. The 
site dilution problem of two-dimensional quantum antifer- 
romagnets has been extensively studied. [1,2,3,4,5] In this 
paper, we treat a problem of the difference between Zn and 
Mg substitutions in La2Cu04. [1] 

As far as the present authors know, on the impurity prob- 
lem of square-lattice Heisenberg antiferromagnets nonmag- 
netic impurities have been treated just as static vacancies 
in spin systems in all previous theoretical studies . However, 
in an experiment on La2(Cui_ K Zn(or Mg) x )04, Cheong et 
al. reported that the initial suppression rate of Neel tem- 
perature 



^ Ris = - lim 



1 dXjsi 
o Tn dx 



(1) 



shows impurity ion dependence, i.e., iijg = 3.4 for Zn 
substitution and -R[Jf s ' ) = 3.0 for Mg substitution. [1] This 
result indicates that the nonmagnetic impurities are not a 
role of static vacancies in the spin system. Cheong et al. also 
measured the orthorhombic-to-tetragonal transition tem- 
perature To^t of La2(Cui_ x Zn(or Mg) K )04, and found a 
faster increase of To^t of Zn substitution than Mg sub- 
stitution. On the basis of this observation, they concluded 
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that Mg is a better substantial ion than Zn for realizing 
spin dilution systems. [1] 

We write the Hamiltonian of the parent compound with- 
out impurities as 



H — J E Si ■ Sj + J± ^ Si ■ Sj, 



(2) 



<«.i>il 



where i and j denote lattice points on a simple cubic lattice, 
(i, are nearest- neighbor pairs of sites on a square-lattice 
plane, and (i,j)±_ are inter-plane nearest neighbor pairs. 
Although inter-plane exchange couplings in La2Cu04 are 
more complex than those in the above model, we simply 
treat the three dimensionality within the form of the above 
Hamiltonian. An "ideal" nonmagnetic impurity is intro- 
duced at site O, then the spin system is perturbed by 



K = J E 

(OJh 



Sj Jx 



(Oj)i 



(3) 



The Hamiltonian of the diluted spin system is given by 
Htot = H + YJq H'q, where the sum runs over Nx lattice 
points chosen at random from the total N lattice points. 
On the basis of the Tyablikov decoupling approximation, 
MacGurn calculated configuration averaged Green's func- 
tions with self-energy corrections within first order in x, and 
estimated Ris for a square lattice case, J± — > 0, and for an 
isotropic simple cubic lattice case, J± — J. [6] MacGurn's 
results are i?is( J± — > 0) = ir and i?is(Jj_ = J) = 1.68. The 
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Fig. 1. Initial suppression rate R\s as a function of interplane coupling 
Jj_. 

former value for the square lattice is close to the experi- 
mental values for La2(Cui_ x Zn(Mg) x )04, but other mech- 
anisms are needed to get a quantitative agreement between 
theory and experiment. One of such mechanisms is a weak 
inter-plane coupling. 

We begin with extending MacGurn's calculation to study 
how Ris depends on J±. The result is shown in Fig. 1. We 
find that Ris decreases monotonically as J± increases. This 
behavior is reasonable because the three dimensionality is 
expected to stabilize antifcrromagnctic ordering. If we as- 
sume J± = 5.0 x 10~ 5 J in La2Cu04, [7] then we obtain 
Ris — 2.67. This value is much smaller than R^'^ = 3.4 
and R ( ™ e) = 3.0, and thus we have to search for an addi- 
tional mechanism which destabilizes the antiferromagnetic 
ordering and makes Ris increase. 

The experimental result of R[s^ ^ R^^ gives us a hint 
in searching for the additional mechanism. A natural way 
to understand this experimental result is that the exchange 
couplings around nonmagnetic impurities are changed de- 
pending on the ion radius of the impurities. The ion radii of 
Zn, Cu, Mg are, respectively, 0.74, 0.73, 0.72 . [1] When a Cu 
ion substitutes for a Zn ion, then the four nearest neighbor 
Cu ions of the impurity moves away from the impurity and 
in the case of Mg they approach the impurity, as shown in 
Fig. 2. (We neglect smaller displacements of other Cu ions.) 
Then, for eight exchange bonds surrounding the impurity 
ion, the linear arrangement of lobes of two d x 2_ y 2 orbitals 
and a p a orbital is changed to a distorted one, which leads 
to the reduction of overlap integral between such two Cu 
sites, i.e., t — ► t — At. The reduction of the overlap integral 
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Fig. 2. Schematic representations of (a) a Cu02 plane, where the 
small closed circles represent positions of oxygens, (b) a CUO2 plane 
with Zn-substitution, and (c) a Cu02 plane with Mg-substitution. 
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Fig. 3. Initial suppression rate Ris as a function of a for 
J ± /J = 5.0 x 10~ 5 . 

also gives J J' = J(l - 2Ai/i) for small At/t. 

In order to incorporate this effect, we write the pertur- 
bation due to an impurity at site O as 

H' ^H' -JaY, Si-Sj, (4) 

{i,j)o 

where a = 2At/t, and the sum runs over eight intra-plane 
nearest-neighbor pairs surrounding the impurity site O. 
Following MacGurn's procedure, we calculate Ris for the 
Hamiltonian H to t = H + J2o H'o- ^ n ^ ne nm it J± - ► 0, we 
obtain Ris — it + Act with a rather large coefficient A = 
32 - 16tt + 5tt 2 /2 ~ 6.4. The calculated result for J± = 
5.0 x 10~ 5 J, which corresponds to La 2 CuC>4, is shown in 
Fig. 3. Using this result and experimental results of -R|g Is ' ) = 
3.0 and i?jg n) = 3.4, we obtain a (Zn ^> = 0.12 and = 
0.06 for La2(Cui_ 2 ,Zn(Mg) a; )04, which means that reduc- 
tions of the overlap integral due to Zn and Mg substitutions 
are, respectively, 6% and 3%. The result of > 
shows that Mg is a better substantial ion than Zn, which is 
consistent with the experimental finding by Cheong et al. 
mentioned before. [1] 

In summary, we conclude that the reduction of intra- 
plane exchange couplings around impurities causes a rather 
fast increase of Ris and gives the origin of the difference 
ini?is between La2(Cui_ a; Zn a ;)04 and L^dCui-^Mg^G^. 
As a future problem, it is also interesting to incorporate 
four-spin ring-exchange terms into our model, [8] which is 
expected to get the estimations of a to decrease more or 
less because of frustration effects. [9] 
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